Microscopic model for Feshbach interacting fermions in an optical lattice with 
arbitrary scattering length and resonance width 
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We numerically study the problem of two fermions in a three dimensional optical lattice interacting 
via a zero-range Feshbach resonance, and display the dispersions of the bound states as a two-particle 
band structure with unique features compared to typical single-particle band structures. We show 
that the exact two-particle solutions of a projected Hamiltonian may be used to define an effective 
two-channel, few-band model for the low energy, low density physics of many fermions at arbitrary 
s-wave scattering length. Our method applies to resonances of any width, and can be adapted to 
multichannel situations or higher- pairing. In strong contrast to usual Hubbard physics, we find 
that pair hopping is significantly altered by strong interactions and the presence of the lattice, and 
the lattice induces multiple molecular bound states. 



The crossover of a system of attractive two- component 
fermions from a condensate of loosely bound Cooper pairs 
to a condensate of tightly bound bosonic molecules has 
a long history [1 , and appears in many contexts, in- 
cluding high-temperature superconductivity [2 and ul- 
tracold atoms [3 . Furthermore, near the crossover such 
a system enters the unitary regime where the scattering 
length is larger than any other length scale in the prob- 
lem. The physics of this regime is relevant to many dif- 
ferent fields, bringing together quantum chromodynam- 
ics, holographic duality, and ultracold quantum gases [4 . 
Theoretical study of the unitary regime is generally dif- 
ficult due to the absence of any small parameter. 

Theoretical analysis becomes even more difficult in a 
lattice, as the center of mass, relative, and internal de- 
grees of freedom become coupled, leading to composite 
particles whose properties depend on their center of mass 
motion [5]. Furthermore, strong interactions require the 
inclusion of a large number of Bloch bands for an accu- 
rate description, and this cannot be handled efficiently by 
modern analytical or numerical many-body techniques. 
In addition to general theoretical interest in how fermions 
pair to form bosons in a discrete lattice setting, the study 
of pairing in lattices is of significant practical importance. 
For example, an accurate, systematically correctable, and 
computationally feasible many-body Hamiltonian is nec- 
essary for calibrating fermionic quantum simulators as 
has been done in the bosonic case [6 . 

In this Letter, we describe a general method to derive 
an effective few-band low-energy Hamiltonian for Fesh- 
bach interacting fermions in a lattice from the numerical 
solution of the two-body problem. We call this Hamil- 
tonian the Fermi Resonance Hamiltonian (FRH). This 
method applies to Feshbach resonances of any width and 
for arbitrary scattering length, and all parameters ap- 
pearing in the effective model can be computed micro- 
scopically from the properties of the two-body solution. 
The difference between the bare model and the FRH is 
sketched in Fig. [l] 

The simplest approach to describing Feshbach inter- 



acting fermions is to replace the interaction with a 
pseudo-potential chosen to reproduce the correct scat- 
tering length. When restricted to a single Bloch band, 
this leads to the popular Hubbard model [T which has 
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FIG. 1: (Color online) Schematic of the FRH transformation. 
(a) In a broad Feshbach resonance, all two-particle scattering 
continua (gray shading) are strongly coupled to bare molecu- 
lar bands (solid lines) . Thus all scattering continua are virtu- 
ally strongly coupled, (b) By correctly dressing the molecular 
bands, one obtains a single scattering continuum (gray) plus 
well-separated dressed molecular bands (green), with much 
simpler couplings. This is our efficient, numerically tractable, 
FRH. 



Exact solution for two particles. The basic concept of 
a two-channel model is for an open channel to describe 
scattering between two atoms and a separate closed chan- 
nel to describe bound pairs. While each channel repre- 
sents a single scattering or bound state in the continuum, 
in the lattice it also acquires a band index. Because of 
an inter-channel coupling, the actual molecule is a super- 
position of bands from both channels. 

To treat this problem, we begin with the nonlinear 
eigenvalue problem developed by Biichler [8] for the 
bound state energy at total quasimomentum K, and i?^, 
the coefficients of the closed channel portion of the wave- 
function. As shown in 0[TO], for the bound states of two 
fermions in an optical lattice interacting via a zero-range 
Feshbach resonance in a two-channel model: 



' Xst •) 



K _ r dq^y^ hnm ( \ Inm^ 
'"sK VM/ 



tK (q)/^]?m(q)- 



(1) 

(2) 
(3) 



Here u is the renormalized detuning between the open 
and closed channels, g is the inter-channel coupling, a 
is the lattice spacing, vq is the volume of the Brillouin 
zone (BZ), and the bars in Eqs. ([2][3| denote quantities 
computed in the absence of an optical lattice. We as- 
sume that spin-spin interactions which change the or- 
bital angular momentum are negligible so that the scat- 
tering is purely 5- wave. The optical lattice is assumed 
to be simple cubic with lattice spacing a and potential 
{x,y,z} / overlaps of the di- 

mensionless coupling between the open and closed chan- 
nels are 



(q) 



/ dxdy [V^nq (x) V^m,K-q {y)^ ^ i^) ^sK. (R) , 



where is the number of unit cells, the V^nq (x) are 
Bloch functions with energies ^nq for particles with mass 
m spanning the open channel and (psK (z) are Bloch 
functions with energies for particles with mass 2m 
in a lattice potential 2V spanning the closed channel. 
We have also defined relative r = x — y and center of 
mass 2R = x + y coordinates, and a (r) is a regular- 
ization of the inter-channel coupling. The sum of the 
noninter acting energies of the open channel is denoted 
^^Yn (q) = ^nq + ^m,K-q and the zcro of energy is 
E^-^^ (0). Here and throughout the rest of this work n 
and m are band indices for the open channel, s and t are 
band indices for the closed channel, q is a single-particle 
quasimomentum, and K is the total quasimomentum. 

While Eqs. ([l])-(|3| apply to resonances of any width, 
we focus on the experimentally relevant limit of a broad 
resonance. Narrow resonances are treated in the sup- 
plementary material . A broad resonance in the 
few-body sense is the limit of effective range much 
smaller than the lattice spacing, <C a, and so we 
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FIG. 2: (Color online) Exact two-particle band structures for 
various as in a strong optical lattice. The bound state energies 
for as /a = —5 (purple boxes), —0.1 (red pluses), 0.1 (green 
crosses) , and 5 (blue asterisks) as a function of the total quasi- 
momentum K along a path connecting the high-symmetry 
points in the irreducible BZ F = (0,0,0), X = (-7r/a,0,0), 
M — (— 7r/a, — 7r/a, 0), R — (— vr/a, — vr/a, — vr/a) for a lattice 
with V/Er — 12. The near-resonant points as/a — ±5 lie 
nearly on top of one another, demonstrating universality. 
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vjE^ oo, as/a = —ng'^/Sa^Ejiu = const., to ob- 
tain {8asER/7ra)x^ (£^k)R^ - = 0, where Er = 
h^Ti'^ /2ma^ is the recoil energy and is the 5- wave scat- 
tering length. How can we then obtain the dispersion re- 
lation for fixed a^, i^, etc.? First, fix the energy eigen- 
value and solve the resulting linear eigenproblem for 
l/ttg. This provides exact eigentuples (£^k,<^s,R^) of 
the nonlinear eigenproblem, though it may not be the 
we seek. Second, fix and use the exact tuple nearest 
this value as initialization for a Newton-Raphson itera- 
tion [11] . This two-stage approach converges to a relative 
accuracy of 0.01% in a few tens of iterations [12 . 

Because the eigenequation Eq. ([T]) is invariant under 
translation by any Br avals lattice vector, its eigenvalues 
can be classified according to the total quasimomentum 
and shown as a two-particle band structure. A complete 
classification of the solutions is given in the supplemen- 
tal material [10]. In Fig. |2] we show only the energy of 
the low-energy bound states with completely even par- 
ity under inversion as a function of K for several as /a 
in a lattice with V/Er = 12. We see the appearance of 
several bound states for a fixed 5- wave scattering length, 
in contrast to the continuum. These additional bound 
states arise from the coupling of quasimomenta modulo 
a reciprocal lattice vector induced by the reduced transla- 
tional symmetry. One salient feature is the emergence of 
universality, which is the independence of the dispersion 
from the sign of when \as/a\ becomes large. For non- 
resonant and negative a^/a, picturing the bound states 
as Fermi pairs with twice the mass and twice the polar- 
izability captures the relative spacings between energy 
levels quite accurately, but generally overestimates the 
effective mass of the bound states. This effective mass 
difference is an indication of the coupling between the 
center of mass and relative motion which leads to impor- 
tant properties of the FRH. 

Fermi Resonance Hamiltonian. A promising route to 



3 



describing Feshbach interacting ultracold gases is by a 
lattice projection of a two-channel model in which the 
closed channel appears explicitly in the Hamiltonian. 
However, for a typical broad resonance such models re- 
quire a large number of both open and closed channel 
bands to solve accurately, and so cannot be treated effi- 
ciently. Because the modern context of this problem in- 
volves extremely low temperatures and densities, we can 
look for an effective model valid in these limits which still 
reproduces the correct physics. This is done by replacing 
the model containing couplings between all open channel 
bands with all closed channel bands with a model describ- 
ing an effective resonance between the lowest open chan- 
nel band with a suitable set of effective closed channel 
bands whose properties are set by the two-body solution 
for low densities. This process is displayed schematically 
m Fig. [1] The purpose of this section is to derive such 
an effective Hamiltonian using our two-particle theory. 

We begin by separating our two-particle Hamiltonian 
using projectors L into the lowest open channel band and 
D = 1 — L into all excited open channel bands and all 
closed channel bands. A similar approach was taken in 
Ref. [13^ for the ID case. An analysis analogous to that 
leading to Eqs. ([l}-([2| gives a nonlinear eigenequation 
for the closed channel components of D|?/^) as 
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where the prime on the sum indicates (m, n) ^ (1,1)- 
Here x differs from x in Eqs. ([iJj-Q in that the summa- 
tion excludes the lowest band. We emphasize that the 
renormalization x includes all bands and so the detuning 
and scattering length used in this projected model are 
those of the full (non-projected) and properly renormal- 
ized two-body problem. We call the eigenstates of this 
projected system dressed molecules. Here we label dis- 
tinct eigenstates of Eq. Q by the parameter a. These so- 
lutions share many features of the full solution presented 
above. However, the divergence of the s-wave scatter- 
ing length for the lowest energy completely even parity 
state occurs near Ej^ = 0, indicating that scattering res- 
onances in the lowest open channel band are generated 
by coupling to this state. 

We now assume that, at low temperatures and low 
densities, two particles which are separated by a distance 
large compared to the effective range of the potential 
will remain in the lowest band to minimize their energy. 
When two particles come together they interact strongly 
and populate many of the excited open channel bands as 
well as the closed channel bands. Because it is rare for 
more than two particles to come together, the particular 
populations of the excited states are fixed by the two- 
particle solution. The dressed molecules encapsulate the 
short distance, high energy physics and couple it to the 
long wavelength, low energy physics of the lowest band 



fermions through the Feshbach coupling. The point of 
connection between the few- and many-body physics is 
the two-particle scattering length (equivalently g and u 
for narrow resonances) , which appears as a parameter in 
the equation Eq. Q defining the dressed molecules. 
The FRH is 
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where a]^ creates a particle with spin a in the low- 
est open channel band Wannier state centered at lattice 
site i, Wi{x.); d] ^ creates a particle in the a^^ dressed 
molecule Wannier state centered at site z, Wi,Q;(x, y); 
h-^ is the number operator for fermions in the lowest 
Bloch band; and is the number operator for the a^^ 
dressed molecule state. The set of dressed molecules M. 
which are included dynamically can be determined on 
energetic and symmetry grounds from the two-particle 
solution. At low energies, only the completely even par- 
ity dressed molecule in the lowest sheet [iO is relevant 
to the set as all others either have vanishing on-site 
couplings from parity considerations or are very far off- 
resonance. In order, the terms in Eq. ^ represent tun- 
neling of atoms in the lowest Bloch band between neigh- 
boring lattice sites i and j; the energy Eq = Ei^^/N^ 
of a fermion in the lowest band with respect to the zero of 
energy; tunneling of the dressed molecular center of mass 
between two lattice sites i and j, not necessarily nearest 
neighbors; detunings of the dressed molecules from the 
lowest band two-particle scattering continuum; and reso- 
nant coupling between the lowest band fermions at sites 
j and k in different internal states and a dressed molecule 
at site i. The FRH is a two-channel resonance model, be- 
tween unpaired fermions in the lowest band, and dressed 
molecules nearby in energy. 

We now describe how to calculate the Hubbard pa- 
rameters appearing in Eq. (|6|. The first term is well- 
known from single-band Hubbard models [14 and we 
do not discuss it here. Due to the fact that the solu- 
tions of the projected nonlinear eigenequation Eq. Q are 
also eigenstates of the total quasimomentum, the second 
and third terms may be written as i^Q, = J d'KE^/vo 
and tfj = -J dKe'^<^'-^^^E^/vo. Because the band 
structure is not separable, 7^ '}2i={x y z} ^^i^ dressed 
molecules can tunnel along directions which are not the 
principal axes of the lattice. This is in stark contrast to 
single-particle tunneling in Bravais lattices which always 
occurs along the principal axes. Thus diagonal hopping 
is a key feature neglected in previous approaches. In 
Fig. ([3| (a) we show that diagonal hopping is often of the 
same order of magnitude as the tunneling of open chan- 
nel fermions in the lowest band. The signs and mag- 
nitudes of the tunnelings and particularly the dressed- 



0.01 



0.001 



^^^^^^^^ 

m 



I / 



^000 



111 




states, one can project out higher bands from x as Wc.s 
done for the lowest band, and then include these bancs 
dynamically in the many-body Hamiltonian with reno' - 
malized couplings. In this way, the energetic domain of 
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FIG. 3: (Color online) Hubbard parameters for the FRH. (a) 
The detunings and tunnelings of the completely even parity 
dressed molecule in the lowest sheet as a function of as /a. The 
detuning is negative for < and positive otherwise. The 
solid black horizontal line is the tunneling of a single open 
channel fermion in the lowest band. The nearest neighbor 
dressed molecular tunneling is nearly two orders of magni- 
tude larger than the open channel tunneling near resonance, 
(b) The effective atom-dressed molecule couplings of the com- 
pletely even parity dressed molecule in the lowest sheet as a 
function of as /a. Schematics of the spatial dependence of the 
various coupling processes are shown in the boxes. 

molecule atom couplings are crucially affected by the 
parities of the dressed molecular Wannier functions. We 
stress that only a full lattice solution can reproduce these 
important properties of the Hubbard parameters; the fre- 
quently used harmonic oscillator approximation will fail 
even qualitatively to do so. 

The remaining Hubbard parameter is the dressed 
molecule-atom coupling, which becomes in the limit of 
a broad resonance g/ Ero.^^'^ oo 



9i-k,k- 
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where the renormalized coupling is g^K = 
En/^-Il^-{dxj^/dE^)-K^ and R,, = R, - R,. 
We emphasize that has only implicit dependence 
on the divergent parametera/£^i?a^/^ through g^-K and 
so remains finite, see Fig. pfb). As ^c^k g/a^^'^, the 
transformation to the FRH has the effect of narrowing 
the resonance. In Fig. |3|b) we also see that the on-site 
coupling ^000,000 is the dominant energy scale of the 
problem for large a^/a, and that off-site couplings can 
also be large compared to other Hubbard parameters 
such as the open channel tunneling. Atoms which do not 
lie along a principal axis can pair to form a molecule, 
but this effect is much weaker than diagonal tunneling 
for the completely even parity dressed molecule. 

In the derivation of the FRH we use only the bound 
states of the projected problem and neglect scattering 
states in higher bands. This captures the scattering 
states in the lowest band and nearby bound states, but 
will fail to capture the physics at higher two-particle rel- 
ative energy where scattering states in higher bands can 
play a role. In order to accommodate these scattering 
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n ot the FRH caii^jje extenclecl art 
more dynamical ^elds. Within the/cpnfines of 
annel model arrcrthe constraint of low ener- 
is an accurate representation of the many- 
[]^JlAn]|lj^nian. However, intrinsic three-body pro- 
cesses which are not captured by the two-channel model 
play a role at higher density and lead to corrections to 
the FRH. A discussion of these three-body processes is 
outside the scope of this paper. 

In summary, we have studied the bound state proper- 
ties of two Feshbach interacting fermions in an optical lat- 
tice at a range of scattering lengths and quasimomenta. 
The bound states of a projected Hamiltonian were used 
to identify a numerically tractable, efficient Hamiltonian 
for a low density many-body collection of lattice fermions 
at arbitrary scattering length and low energies, the Fermi 
Resonance Hamiltonian. Our results provide the appro- 
priate starting point for future investigations of strongly 
interacting lattice fermions. 

We acknowledge useful discussions with J. L. Bohn, 
H. P. Biichler, C. W. Clark, D. E. Schirmer, D. M. Wood, 
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supported by the Alexander von Humboldt Foundation, 
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DERIVATION OF THE NONLINEAR 
EIGENEQUATION. 

Here we review the derivation of Eq. 2 of the main 
text for the bound states of two fermions in an optical 
lattice interacting via a zero-range Feshbach resonance. 
All quantities have their same meaning as in the main 
text. The starting point of our analysis is a two chan- 
nel model with the open channel spanned by states of 
two Fermions in different internal states with equal mass 
m and the closed channel spanned by molecular states 
with twice the Fermionic mass and twice the Fermionic 
polarizability. We describe their interaction via an inter- 
channel coupling g which couples the pair of open chan- 
nel Fermions to a closed channel molecule at the center 
of mass and a detuning u between the two channels. This 
gives rise to the scattering amplitude 



1 



l/ag -\- ik -\- r^k"^ 



(8) 



with s-wave scattering length = — 2/i^^/47rfi^z^ and 
effective range rs = nfi^ / fi'^g'^ . Here fi is the reduced 
mass and k the incident wave vector. 

Denoting the wave function of the two Fermions in 
the open channel as ?/^(x, y) and the wave function of 
the closed channel molecules as ^(z), the two-channel 
Schrodinger equation in position representation is 

[E - i^o(x) - ^o(y)]V^(x,y) = gJdza{r)^{z)S{z - R) , 
[E-uo- H^iz^z) = gjd^dya{r)^{^,y)6{z - R) . 

In this expression ^o(x) = — + ^(x) is the single 

particle Hamiltonian for the open channel and Hq^ (z) = 
— ^V2 + 2V(z) is the single particle Hamiltonian for the 
closed channel. The subscript in denotes that this 
is a bare detuning entering the microscopic theory which 
is related to the physically observable detuning v in the 
limit as the regularization cutoff A ^ oo. Additionally, 
we note that the Feshbach regularization a(r) 5{y) 
in the limit as A oo, where 5{y) is the Dirac delta 
function. 

The open channel solution with total quasimomentum 
K may be parameterized as 



V^K(x,y) 
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where is the total number of unit cells in a 3D lattice 
with periodic boundary conditions and V^nq(x) is a Bloch 
eigenfunction of the single-particle Hamiltonian. As in 
the main text, quantities denoted in bold represent three- 
component vectors, e.g. n = (n^;, n^y, n^). Similarly, we 
parameterize the closed channel wave function as a sum 
over Bloch states computed for twice the mass and twice 
the polarizability ^sk(z) as 



6K(z) = ^i?f0sK(z). 



(10) 



Inserting these expansions into the two-channel 
Schrodinger equation yields 



[E^-v,-EM^]Rf 
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Formally solving the first of the two equations with a 
Green's function and inserting into the second equation 
gives 



^0 



K 



^0 ^mnt £;K-£;S^^(q)+i77 * 



This expression diverges in the limit A ^ oo, as is well 
known for two-channel theories involving a pointlike bo- 
son [2 . We remove this divergence through renormal- 
ization, replacing the bare detuning uq with the physical 
detuning u by subtracting the infinite constant yield- 
ing Eq. 2 of the main text. The divergent parts of Eq. 3 
in the main text cancel and we may safely take the limit 
A ^ oo. 

Following Ref. [1 , we use the regularization 



a(r) 



Jv{A) 



(2^ 



_ikr 
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where the cubical volume v{A) = vqA^ is centered around 
k = with vq the volume of the BZ. We also define a shell 
summation over bands with shell parameter 5, 
as the summation over all band indices n and m less than 
or equal to S with at least one of the band indices being 
S. The correct limiting procedure to obtain x^(K) in the 
limit of an infinite summation over bands and vanishing 
short-distance cutoff is 



lim 

A^oo 



lim xftiEK) 

S^OO 



(12) 



The A limit is taken using the asymptotic scaling relation 

[x^(i?K)](A) = C3t/A + x^(i5K). (13) 

One may be concerned that the scaling relation 
Eq. ([13]) only holds for of the full model and not 
for x^i^ the projected model. To show that this is not 
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the case, we note that in the hmit of an infinite number 
of unit cehs N ^ oo the overlaps h^^{q) may be written 
as products of ID overlaps h^j^i^) form 



r+r'+//27r 



mK-q^M-sK 



(14) 



,2q- K - f + 2Tr(r -r')^ 
xrect(^ 

where / is an integer multiple of 27r which shifts K — q 
into the BZ, rect(x) denotes the rectangle function, and 
the vectors Cnq denote the Fourier expansion of the open 
channel Bloch functions as 



-I'Kirx I a 



(15) 



Similarly, cm;sK denote the Fourier coefficients of the 
closed channel Bloch functions, differs from in 
the exclusion of all terms with n = m = 1. However, 
provided that A is large enough to capture the support 
of the vectors c 



with n 



1, 



Eq. (14) demonstrates 
that these terms are no longer functions of A. Thus, the 
scaling relation Eq. (13) also holds for x^- Similar argu- 
ments show that the same scaling holds for x^ when any 
finite number of open channel bands have been projected 
out. 



CLASSIFICATION OF THE TWO-PARTICLE 
BOUND STATES 

For the simple cubic lattice we consider the Hamilto- 
nian is invariant under refiection in any Cartesian direc- 
tion: H{0^x.) = H(yi') = H(x.) where = ^t^x is related 
to X by changing the sign of all coordinates in some set 1Z: 
Xj —Xj^ j G IZ. Because the generators of refiection 
and translation do not commute we cannot find simulta- 
neous eigenfunctions except at high-symmetry points of 
the BZ. However, the fact that the Hamiltoninan com- 
mutes with both operators implies that parity transfor- 
mations yield relationships between degenerate sets of 
Bloch functions. In particular, for the given lattice po- 
tential, the invariance under the refiection symmetry gen- 
erated by O-jz implies that the Bloch functions transform 
as 



^nq(xO 



n( 



1)"^ + Vn,q'(x), 



(16) 



where we begin indexing the bands from 1. We can 
thus characterize the bands according to whether they 
are even or odd under inversions by the triple p = 
{PxiPy^Pz)^ where = (—1)"^^+-^. This inversion rela- 
tionship implies that the inter-channel overlaps transform 



^sK 



'(q') 



n( 



1) 



/i"K(q), (17) 



and x^ transforms as 

xsT(^k)= n(-i)'^'^'^x^(^K). (18) 

jen 

It can be proven that this transformation leaves the 
eigenvalues invariant, but the eigenvectors transform 
according to 



R. 



K' 



jen 



(19) 



For a total quasimomentum K all of whose components 
consist of either or — 7r/a, this implies that only molec- 
ular bands which transform identically under complete 
inversions mix. Hence, at these exceptional points of the 
BZ, we can unambiguously determine the parity of the 
two-particle state V^kq;(x, y) by the components of its as- 
sociated eigenvector R^. The parity is then chosen to 
depend only on the eigenstate index a by requiring that 
R^ is a smooth function of K. This construction follows 
that of the ID case studied by Kohn 0, which leads to 
maximally localized Wannier functions. 

With this construction, there is still an undefined 
global phase under inversion that we can fix in the follow- 
ing way. The complete two-particle bound state solution 
is 

*Ka(x,y) = -^[^i?f„</.sK(x)f(x-y) (20) 

^??a/lsK(q)V'nq(x)VmK-q(y)l 
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E 



where A/rq; is a normalizing factor and f(x — y) denotes 
a relative wavefunction for the closed channel which has 
characteristic width a/ A. As the theory has already been 
regularized, we may take A ^ oo with the understanding 
that this relative wavefunction has a probability density 
of 1, and forces the closed channel to contribute only at 
the center of mass. Because of the partitioning of Hilbert 
space into open and closed channels, the normalization 
coefficient is 



Ero?!^ 



fn^-x'{E^/En)-B^. (21) 



Here x'{E) is the derivative of x with respect to E. Using 
the transformation properties under O-jz, we find 



*K'a(x,y) =Pa*Ka(x^y')• 



(22) 



Accordingly, we set = ]\j 



This implies that 



LjenP^' 

the dressed molecular Wannier functions transform as 
W,«(x^yO = Pc.Wia(x,y). 

In Fig. [4] we display the bound state energies at K = 
for a lattice of depth V/Er = 12 classified according to 
their parity. The red points correspond to p = (1, 1, 1), 
the green points to p = (1, 1, — 1) et. eye, the blue points 



-15 



-1 



-0.5 




0.5-0.5 1 



FIG. 4: (Color online) Classification of two-particle bound 
states. The bound state energies at K = for a lattice of 
depth V/Er — 12 are classified according to their parity 
and sheet indices. Red corresponds to p = (1,1,1), green 
to p = (1,1,-1), blue to p = (1,-1,-1), and magenta to 
p=( — 1,— 1,— 1). In addition, higher molecular bands for the 
p = (1, 1, 1) level are shown. These give rise to weak scatter- 
ing resonances with the lowest open channel band for as /a > 
and avoided crossings in the higher sheets for as /a < 0. 



FIG. 5: (Color online) Bound state energies for finite width 
resonances. Shown are the lowest energy bound state ener- 
gies for p = (1,1,1) at K = in a strong optical lattice with 
V/Er = 12. The red solid line is rs/a = 0.01, the green 
dashed line is r^/a = 0.1, the blue dotted line is rs/a = 1, 
and the magenta short-dashed line is rs/a = 10. For nar- 
row resonances (large rs) the divergence of as is sharply pro- 
nounced around a narrow energy range, and is shifted down- 
wards from the broad resonance value, compare Fig. [4] 



to p = (1, —1, —1) et. eye, and the magenta points to 
p = (—1,-1,-1). In contrast to the continuum where 
there exists at most one bound state for fixed scattering 
length ttg, there is the possibility of several bound states 
for fixed in the lattice due to the reduced translational 
symmetry. Thus, the parity and the quasimomentum are 
not sufficient to completely describe the states. For a 
fixed 5- wave scattering length a^, we provide two other 
indices which we call the sheet index a G {1,2,. ..,00} 
and the molecular band index s G {1,2,..., 00}. The 
sheet index labels the open channel two-particle scat- 
tering bands, with the convention that the first sheet 
lies below the first band, the second sheet lies between 
the first and second bands, etc. as indicated on the fig- 
ure. The open channel scattering bands are denoted by 
solid stripes. The molecular band index labels eigen- 
states which have the same parity and sheet indices but 
differ in energy. The number of molecular bands ob- 
tained is restricted by the number of closed channel bands 
used to construct x^- Let us define m to be the max- 
imum value of the closed channel band index along any 
Cartesian direction. In Fig. [4j the solid red line cor- 
responds to the completely even parity state computed 
with m = 2 and the red points correspond to the com- 
plete even parity states computed with m = 3. m = 2 
captures the physics well near the lowest open channel 
scattering band. For a^/a > 0, the higher molecular 
bands cross the lowest open channel scattering contin- 
uum at narrow ranges of a^/a, leading to weak scattering 



resonances. For as /a < 0, the higher molecular bands are 
present at extended ranges of a^/a, and avoided crossings 
between these molecular bands can lead to differences 
with lower m computations, see e.g. the third sheet near 
as/a — —0.2. 

FINITE WIDTH RESONANCES 

We now turn our attention briefly to the case where 
g/Ena'^^'^ and v/En are finite. In this case we rearrange 
Eq. 2 of the main text to read 



E 

t 



(^K 



sK)<^st 



R^ = i^Rf (23) 



which is an ordinary eigenvalue equation for the detuning 
u when £'k and g are treated as fixed. We note that g 
cannot be scaled out of this equation as the molecular 
band energies E^ depend only on the lattice strengths 
and masses and not on the resonance width. The solution 
of this equation for various tb and K = is shown in 

Fig.m 

We characterize the width of the resonance in terms of 
the experimentally measurable effective range tb which 
defines the width as g/Ena^^'^ = ^/l6a/7^^rB^ For nar- 
row resonances with large tb only the lowest resonance 
can be seen, and as /a is greater than 1, corresponding to 
strong interactions, only in a very narrow energy range. 
As the resonance becomes broader the energy range over 



8 



which the system is strongly interacting widens, and we 
begin to see resonant behavior near higher scattering 
continua. AdditionaUy, the positions of the narrow res- 
onances are shifted downwards in energy with respect 
to the broad resonances, eventuahy becoming the free 
molecular band energies. We note that the broadest res- 
onance shown is in fact narrower than typical broad res- 
onances found in the experimentally relevant ultracold 
atomic systems, but differs from the infinitely broad res- 
onance results by at most a few percent. This justifies 
our use of the infinitely broad resonance limit in the main 



text. 
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